Quantum phase transition in easy-axis antiferromagnetic integer-spin chains. 
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Antiferromagnetic Heisenberg integer-spin chains are characterized by a spin-liquid ground state 
with no long-range order, due to the relevance of quantum fluctuations. Spin anisotropy, however, 
freezes quantum fluctuations, and the system is magnetized in the presence of a sufficiently large 
easy-axis anisotropy. We numerically investigate the case S = 1, by means of the density-matrix 
renormalization group, and find that the freezing of the spin liquid into a Neel spin solid, with 
increasing easy-axis anisotropy, is a continuous quantum phase transition. Numerical evidence 
indicates that the transition is not in the two-dimensional Ising universality class. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.40.Mg 



Quantum fluctuations play a crucial role in antifer- 
romagnetic (AFM) Heisenberg spin chains, as they are 
strong enough to suppress magnetic long-range order 
at zero temperature. Despite the common absence of 
a spontaneous staggered magnetization in the ground 
state, the physical properties are, however, different for 
integer and half- integer spin 5* [J. Indeed, half-integer- 
spin chains are characterized by a gapless excitation spec- 
trum, power-law spin-spin correlations, and a divergent 
linear response to a staggered magnetic field, and are 
therefore usually considered as nearly ordered. The spec- 
trum of integer-spin chains is instead gapped in the ther- 
modynamic limit, with a finite linear response to stag- 
gered magnetic fields. The ground state is a spin liquid 
with finite spin-spin correlation length, and the gap to 
the first excited state is usually called Haldane gap. 

In this paper we present a numerical investigation of 
integer-spin (namely S = \) chains, for which the pres- 
ence of different energy scales makes a unified theoretical 
approach problematic 0]- This difficulty is mirrored by 
a richer variety of physical regimes with respect to half- 
integer-spin chains. Physical effects which freeze quan- 
tum fluctuations are expected to destroy the spin-liquid 
phase, possibly leading to a magnetization. In a previous 
paper |3| two of us focused on the freezing of the spin liq- 
uid induced by a staggered magnetic field, and found that 
the nature of the spin-liquid ground state entails a pecu- 
liar non-monotonic dependence of the staggered magne- 
tization on the system size in the region of small fields 
(i.e., whose characteristic energy scale is smaller than the 
Haldane gap), and that the full freezing of quantum fluc- 
tuations into a polarized Neel "crystal" is only achieved 
at fields larger than the exchange coupling constant, with 
a large crossover region at intermediate fields. 

Here, as another effect suppressing quantum fluctua- 
tions, we consider an anisotropy in the AFM coupling. 
Since the Haldane gap is small with respect to the ex- 
change coupling, a small (but finite) anisotropy is ex- 



pected to freeze the spin liquid. Thus a zero-temperature 
phase transition to an ordered phase should take place 
upon increasing the anisotropy. We will focus on easy- 
axis anisotropy: When the anisotropy is large, the system 
approaches the Ising limit, magnetized along the easy 
axis. Since the two limiting fixed points, which corre- 
spond to the disordered and fully ordered phases, are 
characterized by different dynamical exponents z, the is- 
sue arises of determining the evolution from the z = 1 
(relativistically invariant) description of the spin liquid 
within the non-linear a model in the isotropic case, to- 
ward the z = Ising limit. 

We base our analysis on the infinite-size density- matrix 
renormalization-group (DMRG) approach Q, which al- 
lows for a very accurate description of the ground-state 
properties of one-dimensional systems. We provide a 
complete description of the freezing of the spin-liquid 
with increasing anisotropy, measuring both non-universal 
features, such as the spontaneous staggered magnetiza- 
tion and the location of the transition point, and univer- 
sal properties, such as the values of the critical exponents 
and hence the universality class of the transition. The 
analysis of an easy-plane anisotropy, which excludes the 
appearance of an order parameter, leaving the way open 
to a Kosterliz-Thouless phase transition, is technically 
and numerically more involved, and is currently under 
investigation 5] . 

Our starting Hamiltonian is a simple modification of 
the Heisenberg (exchange) Hamiltonian to describe a sys- 
tem with easy-axis anisotropy, 
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where Sr — [8^,8^,8^) is the spin operator on site r, 
8 is integer, and L is the number of sites in the chain. 
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J > is the AFM Heisenberg coijpling, < p < 1 is the 
easy-axis anisotropy parameter [Q], and H is the amph- 
tude of the staggered magnetic field along the z (easy) 
axis, which is coupled with the staggered magnetization. 
As it is customary and convenient in DMRG calculations, 
we assume open boundary conditions. The two limiting 
cases p = \ and p = Q correspond to the Haldane spin 
liquid and Ising antiferromagnet respectively. It is worth 
noting that the model is invariant for p —p, so 
that our results equally apply to a system with a ferro- 
magnetic coupling on the xy plane. 

As discussed in we find it more convenient to per- 
form a transformation on the Hamiltonian which al- 
lows us to make use of the standard implementation of 
DMRG 0. After a local rotation of the reference frame, 
S'^ = i-iyS^, = i-iyS^, Sy = Sy, the Hamiltonian 
is re-cast in the "ferromagnetic" form 

L-l 
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hY^S^, (2) 
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where 5± = ± iSy are the rising and lowering oper- 
ators, and the external field H is now coupled to the 
uniform magnetization 5^ — J2r^r- The easy-axis 
anisotropy corresponds to \p\ < 1, whereas the case 
\p\ > 1 corresponds to an easy-plane anisotropy. In the 
following we refer to the properties of the transformed 
model described by the Hamiltonian and in partic- 
ular to the uniform magnetization along the easy axis, 
(5^). The correspondence with the original AFM model 
is straightforward. 
Differently from the case discussed in Ref. (J], where 
the magnetization was forced by the external field, a 
spontaneous magnetization mo is expected to appear 
now, as quantum fluctuations are frozen with increasing 
easy-axis anisotropy. Indeed, the limiting results, mo = 1 
for p = (Ising limit), and mo = for p = 1 (isotropic 
limit), suggest that mo should decrease with increasing p. 
However, contrary to the case of half-integer spin, where 
the system is critical at p = 1, and is therefore sponta- 
neously magnetized for arbitrarily small anisotropy, here 
we expect that the system persists in the magnetized 
phase only up to a critical value p = Pc < 1, and the 
spin-liquid phase extends over a finite region Pc ^ P ^ 1- 
Since, as we discuss below, the transition is continuous 
(of the second order) , the issue arises whether the proper 
dimension of the effective model at the phase transition 
equals or not the effective dimension of the relativisti- 
cally invariant spin- liquid phase (de// = d-^z = 2, where 
d = 1 is the spatial dimension and 2: = 1 is the dynam- 
ical exponent). It is worth noting that the two phases 
(magnetized and spin liquid) are both characterized by a 
gapped excitation spectrum, although of different origin. 



The gap in the easy-axis phase is related to the discrete 
symmetry, whereas the Haldane gap in the spin-liquid 
phase is of a more subtle nature. 

To characterize the ordered phase and the transi- 
tion, we compute the magnetization M(J, p, H; L, Nk) = 
{[S^]'^), where (•) denotes the expectation value over the 
DMRG approximate ground-state wave function, for the 
system of linear size L, with A'^^ states retained in the 
density-matrix truncation procedure, and the superscript 
T indicates the truncated operators defined in the ap- 
proximate basis. We consider the magnetization per site 
in units of S*, m = M/{SL) (we always take S = 1), such 
that < m < 1 for H > 0. m depends on J and H only 
through h = H/ J. Starting from m(p, h\L, Nk), we per- 
form three subsequent extrapolations [7| to extract the 
spontaneous magnetization and the linear susceptibility: 

(i) We extrapolate the finite-size results to the thermo- 
dynamic limit L ^ 00 for any given number of states Nk 
and field h by fitting the large-L behavior of the mag- 
netization as m(p, h; L, Nk) = 'rn{p, h; Nk) + A/ L, where 
the long-living corrections are 0{1/L) due to the open 
boundary conditions. According to our previous results 
13, the linear fit in 1/L at small field is accurate only 
if the size L of the chain is larger than roughly twice 
the characteristic length scale ^ of the spin-liquid state 
(^ = 6.03 in units of the lattice spacing, for S' = 1 Q). 

(ii) Once m(p, h] Nk) is obtained, we carefully check the 
convergence with respect to the Hilbert space, by taking 
the extrapolation to the limit Nk ^ 00 in the region 
where m displays a clear asymptotic dependence on 1 /Nk 
0, m{p,h-Nk) = m{p,h) +B/Nk. Up to Nk = 100 
states were needed for the smallest field h = 0.0001, to 
obtain the asymptotic behavior. We emphasize that this 
procedure represents a direct measure of the convergence 
to the limit of infinite Hilbert space of the expectation 
values of the different observables, and is not equivalent 
to fix a small truncation error for the density matrix. 

(iii) We finally obtain the spontaneous magnetization mo 
and the linear susceptibility x{p) by means of the extrap- 
olation for small magnetic fields m(p, h) = mo{p)+xip)lT'. 

All the above extrapolations are carefully checked for 
convergence, e.g., with respect to various extrapolation 
schemes, and the least-square linear fits used thereby are 
chi-square tested. The numerical values that we report 
below are affected by an error on the last digit. The error 
bars for the data points reported in the various figures 
are always smaller than the symbol size. 

The change in the ground-state properties with de- 
creasing p is already witnessed by the behavior of 
m{p, h; L, Nk) at fixed Nk and ft, as a function of the sys- 
tem size L (not shown). Due to the presence of a finite 
characteristic length, a non-monotonic behavior charac- 
terizes the spin-liquid phase sd p = 1, with a bump for 
L ~ ^ 3] . Reducing p gradually leads to a monotonic be- 
havior as a function of L, with much larger bulk {L — > 00) 
magnetization. The study of the spontaneous magneti- 
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FIG. 1: The spontaneous magnetization m as a function of 
the anisotropy p. The different data points refer to different 
values of the lowest magnetic field hmin used in the small-field 
extrapolation: hmin = 0.001 (black dots), 0.0005 (full trian- 
gles), 0.0001 (open circles) , and 0.00005 (open triangles) . The 
inset displays the same quantity in log-log plot to emphasize 
the power-law behavior. The dashed line is the fitting curve 
mo = a {pc — p)^ (see text). 



zation, reported below, confirms that the change of be- 
havior is associated with the change in the ground state 
from spin liquid to spontaneously magnetized. 

The magnetization m(p, h) increases with decreasing 
p at each fixed value of h, and naturally increases with 
the magnetic field h at each p. The dependence on h 
is weak at p ~ 1 and vanishes as p 0. The critical 
behavior is evident in Fig. ^ where we report mg {p) as a 
function of p, obtained through a zero-field extrapolation 
^|. Except for a narrow region around the critical point 
10l |. the data are well described by mo = a (pc — p)^ , for 
p < Pc, and mo = for p > Pc, with a = 1.09, /? = 0.139, 
and Pc = 0.845, as obtained by numerically minimizing, 
with respect to pc, the least-square fit of log mo, which 
yields an "optimum" value of a and (3 for any given pc- 
The inset of Fig. ^shows that 11 extrapolated points fall 
on the fitting curve. 

We can also extrapolate the linear susceptibility x{p) 
for /i — > (see Fig. |21), finding the power-law behavior 
X^^ = b{pc - pY' for p < Pc, and = b'{p - PcV' for 
p > Pc, once again except in a small region close to the 
critical point. Using for pc the value obtained from the 
spontaneous magnetization, pc = 0.845, and performing 
a least-square fit of logx~^, we find b = 41.1, 7 = 1.71, 
b' — 1.30, and 7' — 1.70. The numerical agreement of 
the two indices 7 and 7' within numerical accuracy is 
a check of the reliability of our results. When studied 
over the entire parameter range, as a function of p, xip) 
obviously vanishes in the Ising limit {p ^ 0), where the 
magnetization saturates (mo — > 1), and has a finite value 
in the Haldane spin- liquid phase at p = I, which increases 



FIG. 2: Inverse susceptibility a function of p in a log- 

log scale. Notations are the same as in Fig. 1. The dashed 
lines are the fitting curves = b{pc — p)'^ for p < pc, and 
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for p > Pc (see text). 
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FIG. 3: The magnetization at the critical point p — pc as 
a function of the magnetic field. The inset shows the same 
quantity in log-log scale to emphasize the critical behavior. 
The dashed line is the fitting curve m{h,p = pc) = ch^^^ (see 
text). 



with increasing S" (x = 18.50, for = 1 [g). 

The exponent 5 which characterizes the behavior of 
the magnetization as a function of the external field h 
at criticality, m{h,p = Pc) = ch^^^ , can be calculated by 
means of the scahng law S — 1 + j/f3. Substituting the 
values of 7 and (3 found above we obtain S = 13.4. We 
can compare this value with a direct measure of S. Fig. 
I^lshows the data for m{p = Pc, h), with pc = 0.845, as a 
function of h. The log-log plot shows a linear behavior 
over three decades with a least-square exponent 5 = 13.0 
in agreement with the one obtained from the scaling law 
within numerical accuracy, and c = 0.938. 
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FIG. 4: Correlation length ^ as a function of the anisotropy 
p. The dashed line is a guide to the eye. 



Ising universality class (although quite close to it) . Thus 
in the evolution from the spin liquid (where de// = 2) 
to the Ising antiferromagnet (where d^ff = 1), a fixed 
point is reached with a different effective dimensionality 
deff — 2.3 (i.e., a dynamical exponent z = 1.3). This 
result indicates that the lattice-renormalization effects 
and/or the less relevant coupling constants neglected in 
Ref. ^1 , play a role, even if quantitatively small, so that 
the critical exponents are close to, but different from, 
those of the two-dimensional Ising model. 
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lani. We also acknowledge financial support of INFM/G 
through PA-GO-4 and of MIUR Cofin 2001, prot. 
2001023848. 



Finally, by calculating the spin-spin correlation func- 
tion along the easy axis Cif(r) = {S%^^^_^^[Sl^^r^^2rV) 
(r = 0,1,2,...) [ll|, we determine the exponents v,i'' 
which characterize the divergence of the correlation 
length ^, given by ^ — g{pc — p)~^ for p < pc, and 
^ = g'{p — Pc)~'^ for p > pc- The correlation length 
is obtained by fitting the long-distance behavior of the 
correlation function as Cif(r) = Ae~^^^ + B, where 
B = Cxir — > oo). ^ (see Fig. 0J. By fitting the data 
for log^ as a function of log |p — Pc|, assuming pc = 0.845 
as determined from the fit of the spontaneous magneti- 
zation, we find i' = 0.86, g = 0.45 (fitted with 11 points) 
and J/' = 0.87, g' = 1.04 (fitted with 12 points). It is 
worth noticing that the correlation length ^ as a func- 
tion of the anisotropy parameter p is rapidly reduced to 
a quantity of the order of the lattice spacing as p is re- 
duced on the Ising (ordered) side, and reaches the finite 
value ^ ~ 6 in the Haldane spin- liquid state, at p = 1. 

Using the scaling laws, we can determine the effective 
dimension at the phase transition c?e// = (2/3 + 7)/j^ = 
2.3. We have therefore numerical evidence that the fixed- 
point effective model is not relativistically invariant. In 
the flow from the spin-liquid to the fixed point, along 
with the various critical exponents, also the dynamical 
exponent z is therefore corrected in a non-trivial way. 

In summary, we have studied the role of easy-axis 
anisotropy in the S = 1 AFM spin chain by means 
of DMRG. Our calculations show that the freezing of 
the quantum fluctuations due to increasing easy-axis 
anisotropy at small magnetic field is much alike the one 
obtained with increasing field In particular the non- 
monotonic dependence of the magnetization m on the 
size of the system L observed for the isotropic chain is 
gradually transformed into a monotonic increase. The 
extrapolation to zero magnetic field reveals a quantum 
phase transition from a spin liquid to a spin solid at a 
finite value of the anisotropy parameter, pc — 0.845. The 
values of the exponents indicate that the phase transition 
falls within a different universality class than the d = 2 
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